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ABSTRACT 


In this thesis, two important dynamic aspects of tethered 
satellite system are investigated. In the first part, attention 
is focussed on payload raising potential of dumb-bell tethered 
satellite configuration executing pitching llbration in circular 
as well as elliptic orbits. The effect of librations on the 
maximized gain in orbit size subsequent to tether cut-off at 
optimal location is examined. The study clearly demonstrates the 
usefulness of having the dumb-bell tethered system in highly 
elliptic orbit for substantially enhancing the tether’s payload 
raising capability. In some extreme cases, e.g. with e = 0.5, the 
increase in the payload major axis can go up to as high as ~ 50 

times the length of the tether, much more than the corresponding 
figure reported in the recent literature. 

In the next and final phase, the ln-plane transverse 
vibrations of the two tethers connecting the three-body 

constellation in its stable equilibrium configuration are 
investigated. The natural frequencies of vibration in the first 
few Important modes are obtained using finite element approach. 
The results obtained here may be of considerable general 
importance for a preliminary design of various components or 
subsystems in the tethered constellation. Besides, the analysis 
presented may be found useful for extracting information on 
tension in the vibrating tethers connecting the satellite 
constellation. 



CHAPTER I 


INTRODUCTION 


1.1 PRELIMINARY REMARKS 

Advent of the space shuttle and the space station have 
presented a vide range of possibilities for space exploration and 
exploitation. One approach to this end is the concept of Tethered 
Satellite System (TSS). However, the idea is not new. Indeed in 
189S, Tsiolkovsky suggested linking large masses by a long thin 
string to exploit weak gravity gradient forces * _ 

The possibility of using tethers In space has generated 
considerable interest. Tethered Satellite Systems are general 
candidates for enhancing future shuttle capabilities. They allow 
the deployment and retrieval of payloads down to orbits which 
cannot otherwise be considered due to high drag decay. 

Operations with tethered vehicles were successfully carried 

2 

out during two Gemini missions using a short tether. The manned 
spaced vehicles Gemini XI and Gemini XII were tethered to unmanned 
Agena vehicles. The former used a rotating configuration that was 
maintained by the Gemini thruster reaction control system, whereas 
Gemini XII flight test demonstrated a gravity gradient specialized 
configuration. 

This very year l.e., in 1990, it was proposed to fly the 
first tether deployer/retri'eval system in space. This mission, 
delayed due to temporary setback to the shuttle programme 



following the challenger disaster, was to fly the Tethered 
Satellite System ’ developed and funded by NASA and the Italian 
National Space Plan (CNR). The first mission is to deploy a 500 
kg subsatellite upward from the shuttle on a 2 mm thick Kevlar 
tether 20 km long, where the gravity gradient forces will 
stabilize the two masses along a vertical configuration. It is 
proposed to investigate the interaction of the tether, having a 
conducting core, with Earth’s magnetic and ionospheric 
environment. A second mission is to be flown about two years 
later, using the same deployer system, to tether a subsatellite 
100 km downward (the shuttle will be at 220 km altitude) to 
investigate Earth's upper atmosphere. Several other tether or 
tether related missions are in the planning and development stage. 
The 1990’s should see at least half a dozen such missions which 
will demonstrate specific tether capabilities and increase our 
understanding of their benefits and behaviour in space^. 

Several possible applications of TSS have been proposed for a 
multitude of uses. A few of them are summarized below: 

(1) Low altitude scientific experiments : A subsatellite may 
be deployed into the upper atmosphere from the shuttle to provide 
scientists with a means of measurements of the physical properties 
of atmosphere. It would also permit long term observations of 
various phenomena in the lower thermosphere. 

(2) Deployment of payloads into new orbits or retrieval of 
satellites for servicing . 

(3) Antenna for radio astronomy and low frequency 


communications . 



(4) Power Generator : A long Insulated conducting tether 
moving at a high speed cutting Earth's magnetic lines of force 
induces e.m.f. Such a system could be used as an emergency power 
generator. 

(5) Micro gravity experiments : A tether connected system 
aligned with the local vertical can produce an artificial gravity 
for both the end bodies. The tension in the tether caused by the 
gravity gradient pulls these bodies towards their system centre of 
mass (CM). The weak gravity level thus produced (0.01 to O.lg) 
makes this system ideal for low gravity experiments. 

(6) Transportation and Space Constellation : A tether can be 
used for cargo transfers or deposition of hazardous facilities to 
a safe distance or even for capturing an artificial comet. 

Besides, there are many other possible applications proposed 
but not discussed here for the sake of brevity. 

The conceptual use of tether in affecting an overall 
propulsion economy for "enlarging” a payload orbit is now well 
established. This is accomplished by launching the space 
transportation system along with its tethered but undeployed 
subsatellite into an orbit. The subsequent deployment of small 
payload through tether followed by its release raises the small 
payload orbit while the shuttle deploying it gets correspondingly 
lowered in this "momentum exchange" manoeuvre. In recent years, 
several Investigators have shown that release and capture of 
payloads through tethers are practical and fully stable 
operations. The underlying concept is as follows: the 
subsatellite when held by a tether deployed upwards has the same 
angular velocity as the main orblter, which implies a subsatellite 
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velocity greater than that needed for the orbit at its altitude. 
Thus if now released by cutting the tether, it will lie on an 
®Hlptic orbit around the Earth having a greater semi-major axis 
than that of the tethered two body system. More explicitly-, the 
three steps involved in this operation are as follows: 

(i) Deployment of the payload to a planned release altitude 
using tether. 

(li) Disconnection of the payload from the tether enabling it 
to rise it to a higher orbit. 

(ill) Reeling back the tether into the shuttle, an operation 
called retrieval. 

First part of the thesis deals with this application. 

Dynamics of the tethered system is rather complicated. It 
Involves three Important aspects of motion. The orbital motion 
refers to the translational motion of the CM of the tethered 
satellite system while the attitude dynamics is concerned with the 
"rigid body" angular motion of the system about the CM. Unlike the 
rigid body system, here we also have a flexible tether which is 
responsible for sharing the energy in its longitudinal as well as 
transverse modes of vibrations. For complete understanding of the 
system dynamic behaviour, in general, all the three aspects must 
be examined simultaneously. This problem gets further complicated 
by the presence of environmental disturbing effects such as 
atmospheric drag, aspherlcity of the Earth, Earth’s magnetic 
field, gravitational pull of other celestial bodies, etc. 

The dynamics of tethered satellite systems has received a 
great deal of attention during the last decade and at present 
fairly complete dynamical models are available. But with a few 
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exceptions, most investigations usually deal with the rotational 
motion of the tether. However, of the various aspects of motion 
of the tethered satellite systems, the associated elastic 
oscillations of the tether have received relatively little 
attention. Later part of the work is devoted to the transverse 
vibrations of the tether in the three body system connected by a 
pair of tethers. 

1.2 A BRIEF REVIEW OF THE LITERATURE 

A fairly large number of models used by various investigators 

in an attempt to obtain some insight into the complex dynamics of 

6 7 

the system have been summarized by Mishra and Modi ’ . Here, we 

present literature review separately for the two topics: release 
of the payload of a two body system and the transverse vibrations 
of the three body tethered constellations. 

1.2.1 Release at Payloads 

The general concept of releasing a tethered payload from a 
shuttle is shown in Fig. 1. Uhen deployed upwards the payload 
travels faster than the circular orbit speed required at its 
altitude. On the other hand the shuttle possesses a speed lower 
than what is needed to sustain its orbit at this altitude without 
the tether. If the tether is now severed, the payload would be in 
a transfer orbit to a higher apogee while the shuttle would drift 
to a lower perigee. 

A 

Von Tiesenhausen reported that in order to have the released 
payload ending up in a circular orbit, the initial shuttle orbit 
needs to be elliptical. 

Bekey 9 reported that for a two body tethered system in a 
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circular orbit the released payload would follow a different path 
with its perigee at its release point and the apogee on the 
diametrically opposite side of the Earth. If the release occurs 
while the payload is librationally stable and along the local 
vertical, the altitude of the apogee would be higher than the 
perigee ~ seven times the length of the tether. A suitably 
pitching payload can rise to much higher with the corresponding 
gain now going up to fifteen times the length of the tether. 

Zine Eddine Amler 1 ^* considered a maximization procedure to 
determine the optimum time of release of the payload so that the 
altitude gain for the two body system in a circular orbit on 
tether cut-off Is maximized. He concludes that the gain In apogee 
altitude for the librating tethered satellite system can be up to 
even as high as sixteen times the tether-length. 

1.2.2 Space Constellations 

Until recently, the Interest of the investigators in tethered 
system applications was focussed on two body tethered systems. 
There are, however. Interesting situations involving three 
orbiting bodies connected by tether(s), such as the space station 
based tethered constellations and transportation of a cargo from 
one end body to the other. A space constellation is defined as a 
fixed arrangement of three or more masses connected by tethers. 

Since the concept is relatively new, there are few dynamical 
studies in this area. The first dynamical formulation of a three 
body tethered system was attempted by Liu 11 . He considered the 
attitude dynamics in the orbital plane during the mass transport 
from one end body to the other. His formulation led to six 
non-linear ordinary differential equations in the six variables. 
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Since there can be only three independent generalized coordinates 
in this case, three constraints were also obtained. In an 
alternate formulation, he utilized four generalized coordinates 
with just one constraint. Unfortunately, in either case the 
author found it difficult to solve, even numerically, the 
complicated system of equations subjected to the algebraic 
constraints . 

A more efficient formulation for the three body tethered 

12 

systems was presented by Lorenzini et.al. Neglecting the mass 
of the tethers the frequency of oscillation of the line joining 
the end bodies with respect to the system centre of mass was found 
to be approximately -1)0, where R£ is the ratio of 
mass of the middle body to the total mass and O is the overall 
system orbital frequency. 

Zlne-Eddine Amier*^ has derived the equations governing the 
two-dimensional dynamics of three body tethered constellations in 
an Earth Orbit. He has also determined the equilibrium 
configurations and stability of motion around the equilibrium 
configurations . 

1.3 PURPOSE AND SCOPE OF THE INVESTIGATION 

From the literature review, it is clear that the attitude 
dynamics, i.e. the swing motion, can play a major role in the 
release of payload and signlf leant ly alter the orbital motion of 
the released payload. On the other hand, the dynamical analysis, 
especially the vibrations of the tethers connecting three body 
constellation is at its infancy. Here, we consider these two 
topics as our main objective for the investigation. First, the 
effects of the various parameters of the two body tethered system 
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on the maximum gain in the major axis of the payload, starting 
from an elliptical orbit, are examined. Next undertaken is an 
analysis of the transverse tether vibrations of three body 
constellations in a vertical stable configuration. 



CHAPTER II 


TWO-BODY TETHER SYSTEM 


2.1 INTRODUCTION 

A subsatellite deployed from an artificial earth satellite, 
and attached by a tether constitutes a two-body tether system. 
During deployment or retrieval, the system centre of mass 
continues to move along a Keplerian orbit. In general, the 
dynamic study covering various aspects of motion is fairly 
involved. However, several of these can be investigated with 
sufficient accuracy by simplified tether models. In our 
formulation, the mass of the tether is taken into account while 
obtaining the total kinetic and potential energies of the tethered 
system. The kinetic energy is calculated as a sum of that 
associated with a translational motion as well as that for the 
relative altitude motion about the mass centre. Thus the dynamics 
is separated in the orbital motion of the centre of mass and the 
relative motion of end masses with respect to the local frame 
connected to the centre of mass ([attitude motion). No 
restrictions on length of the tether or on mass ratio of the 
end-masses are imposed. 
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The system equations of motion are developed using Lagrangian 
formulation assuming the external forces to be only gravitational. 
The equations are solved with a view to study the variation of the 
total energy of the subsatellite as the two-body tether system 
orbits around the earth. This enables us to identify the Instant 
at which the total energy of the subsatellite is a maximum. It is 
at this point when cutt lng-of-the tether maximizes the momentum 
transfer from the shuttle to the subsatellite. Subsequent growth 
in the satellite orbit is assessed in terms of the net Increase in 
the semi-major axis of the payload achieved over that of the 
tethered system. For simplicity, we only consider the pitching 
librations. The out-of-plane librations and tether vibrations are 
ignored . 


2.1.1 Approach to the Problem 

To find out the optimum time of release to maximize the 


"altitude gain", various methods are available in the literature. 
A few of them are: (1) Maximizing velocity of the subsatellite at 
the time of cutting the tether, (2) Finding pitch angle (o) and 
its derivative at every position along orbit and making 
simultaneously ot to be maximum and «' to be zero at the point of 
cutting the tether. All these methods seem to work fine as long 
as Initial orbit of the system c.m. is circular. For elliptic 


orbits, since system equations 


motion 


non-linear 


differential equations, these methods tend to diverge 


eccentricity of the orbit increases. 

In the present Investigation, a new approach has been 
followed. Instead of trying to maximize velocity or pitch angle 
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(«) at the time of cutting, what has been tried is to maximize the 
total energy. As is known, the total energy per unit mass for the 
subsatellite is 

« = -p/2a, 
where, 

N is Earth’s gravitational constant 
a is semi-major axis 

The orbit apogee, can be calculated from the instantaneous 
position and velocity vectors of the subsatellite at the cut-off 
point. The position and velocity vectors become known once we 
solve the system equations of motion. These equations are 
non-linear and non-autonomous . An approximate analytical as well 
as an "exact" numerical solution has been attempted. 

The approach followed here turns out to be a simple and 
elegant one. It works for both circular and elliptic starting 
orbits and no divergence has been observed. 

2.2 DEVELOPMENT OF THE SYSTEM EQUATIONS OF MOTION 
2. 2. 1 Description of the System 

The system under consideration is shown in Fig. 2. Hasses 
and represent the orbiter (shuttle or space station) and the 
payload (satellite) respectively. The long tether joining and 

H 2 has a mass p t per unit length and length l t . The instantaneous 
centre of mass (CH) of the system is located at a radial distance 
R along the local vertical from the centre of the Earth and its 
true-anomaly is d. Ue consider three sets of coordinate axes. 
The global coordinate system, XYZ, is an Inertial reference frame 
fixed at the centre of the Earth. The other two coordinate 



12 


systems, xyz and x o y o z c &re rotating frames having their origin at 

the system centre of mass (CM). In x y z coordinate system, z 

axis is along the local horizontal, y^-axis is along the local 

vertical and x Q -axis is normal to the orbital plane. In xyz 

coordinate system, however, y-axis lies along the tether line, 

z-axie is normal to y-axis and In the orbital plane while x-axis 

Is normal to the orbital plane. Ue limit ourselves to only 

ln-plane librations. Thus the orientation of the tethered system 

relative to ' x .^ y q ^ q frame is described by pitch angle ot only. The 

pitch angle oi is the satellite rotation about orbit normal 

x -axis. The unit vectors along the three set of axes are 
o 

respectively I, J, K; 1 q , j Q , and i, j, k. The radius vectors 
r^ and ^ denote the position vectors of masses M^ and M 2 from the 
CM of the system. The position vector of any arbitrary point on 
tether, with respect to the CM, is denoted by and can be 

4 4 

determined uniquely in terms of r 1 and • 

Orbital dynamics and attitude dynamics have been considered 

separately assuming the perturbation of the orbital motion due to 

1 2 

attitude motion to be negligible . However, the attitude motion 
of the entire system may have non-negligible Influence on the 
orbit of the payload after release. The other assumptions are: 

(i) Rigid tether, ignoring its vibrations - this part is 
separately considered for space constellations in Chapter 
III. 

(il) The orblter and payload regarded as point-masses . 

(iii) Only external force acting on the system being the Earth's 

the other forces like atmospheric 


gravitational force. 
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drag, effect of Earth’s electromagnetic field, solar 
radiation pressure, gravity forces due to other space 
bodies etc. are ignored. 

These assumptions considerably simplify the mathematical 
model, but the analysis is primarily directed to bring out the 
basic features so essential to study an advance model. 


Z* Z» Z Kinenaiicfi of the System 

Referring to Fig. 2, since CM is the centre of mass, we have 
the following relation. 


W 1 *1 + M 2 ^2 + J ^3 d ® = 0 ( 2 . 2 . 1 ) 

M 3 

And, from geometric considerations, 

*2 ' = *t (2.2.2) 

where |^ t | is 1^, the length of tether. 

If we resolve the vectors in their components along the x, y, 
z axes introducing a spatial variable y^ measured along the 
tetherline from the mass , it allows us to write 

r 3 = (y t - j 

A 

? 2 = r 2 j 

*1 * - r l ! 

Substituting these into (2.2.1), we have 

f -( Bl * fl 2 ) rj * n 2 l t ♦ p t S (y t - r l ) dy t ] J = 

0 


0 
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which yields the following equation 

-CHj ♦ h 2 * m 3 ) r 3 * cm 2 ♦ \ m 3 j i t 


= 0 


i . e . 


H 2 + I H 3 
H 


(2.2.3) 


where, denotes the mass of the tether. 


Similarly, 


M 1 + 7 M 3 
H 


(2.2.4) 


where, the total system mass M = + M^. 

2. 2. 3 Kinetic and Potential Energies 

The total kinetic energy and the total potential energy, each 
can be separated into two parts, one for the orbital motion and 
the other for the attitude motion. For example, 

T = T . + T . 

orb nonorb 

V = V + V 

orb nonorb 

Ue will consider each part separately. For the orbital 
motion, using polar coordinates and & for the CM, we have 

T «,rb ' I" I *.* * *o ® 2 1 <2 - 2 ' 5 > 

— V orb - - ,S_ (2.2.0 

For non-orbital part, we first obtain the energy expressions 
relative to the xyz coordinate system and then use the following 


trans formations 
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r 

O O 

< <*n 

► — 

10 o 

0 cos OL -sin a 

-1 

^ 

V-M > H* > 



0 

i 


0 sin ot cos a 


A 

k 

h. j 


( 2 . 2 . 7 ) 


r " 

I 

J 

»■ sc 

1 0 

0 cos 0 

0 

-sin & 


i 

o 

j K 

K 

L. ^ U 


0 sin 6 

cos d 

; 

k 

o 

b j 


( 2 . 2 . 8 ) 


Also, noting that oi, the angular velocity of xyz frame has 
the following components 


a = (©+«);<*> = = 0 

x y z 


so that 


4 4 -4 

" x r i» r 2 


4 4 *4 

= « X r 2 , r 3 


-4 

= x r. 


1 ‘ 1 ’ ‘2 ~ ~ ‘ 2 ’ ‘ 3 ‘ 3 ’ 

Using these relations to transform the energy expressions, 

after going through fairly straightforward but lengthy algebraic 

manipulations, one can show that 


T = -i J (£ + «) 2 

nonorb 2 


V k = 7 ^ I 1 - 3 co ® 2 

nonorb 2^3 

o 


( 2 . 2 . 9 ) 

( 2 . 2 . 10 ) 


where J represents the moment of inertia of the system around an 
axis passing through CM and perpendicular to the tether and can be 
expressed as 


J = m 1 


( 2 . 2 . 11 ) 



l€ 


where , 


o> 


■ "i ( 


K 2 + I «3 2 


i r 


) ♦«.( 


M 


1 + 1 


m. 


+ w. 


f 1 _ M 2 + I M 3 
1 1 H 




i r 


+ 1 M - 


]} 


In the special case, if one considers to be very large compared 
to M- and H., a reduces to 


a 


m 2 + y n. 


( 2 . 2 . 12 ) 


2. 2. 4, Equations of Motion Before Release 

The equations (2.2.5), (2.2.6) and (2.2.9), (2.2.10) give the 
expressions for the total kinetic energy and potential energy. 
Using the lagrangian formulation, one can show that the equation 

7 

governing the pitch oscillations is : 


or 


_ f 2 e sin g_ 1 

l 1+e cos & J ( } 


1 + e cos & 


cos a sin a = 0 


(2.2.13) 


2.3 ORBITAL MOTION AFTER RELEASE 

In the previous section equations of motion were developed 
for the tether system before release considering both orbital and 
attitude dynamics. However, after the release of the tether we 
have two unconnected point masses, so we need to consider simply 
the orbital dynamics for each of the two separated masses. Since 
our interest lies in the motion of subsatellite, let us denote its 
final orbit apogee and eccentricity by a 2 and e 2 respectively. To 
calculate these we need to know its position and velocity vectors, 
and respectively, at the time of release. 

2 
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2.3.1 Position and Velocity Vectors Just Before Pel ease 
From Fig. 2, we have the vector relation 

*2 = K + *2 (2.3.1) 

where , 


r 2 f 2 1 (2.3.2) 

and the distance r 2 is a function of masses H^, M 2 and as given 

by Eq. (2.2.4), while R^ is obtained from the initial orbit 

parameters a^ and e Q using the following well known relation for 

ellipse in polar coordinates (R q , ©). 

1/R = 1 + e cos & (2.3.3) 

o 

where the semi-latus rectum, 1, is given by 

1 = a (1 - e 2 ) (2.3.4) 

Therefore, if & be the true anomaly at the point of release in the 
original orbit, we have 



✓ 2 . 

“c < l ‘ e o 3 

1 + e cos & 
o 


(2.3.5) 


Also, for the central force motion 

R ^ & - const. = h (2.3.6) 

o 

where, h is specific angular momentum and is given by. 


. 2 . 
h = ^ - 1 ; 

These led to the final results as given below: 


R ~ j- (e sin 6) 
o ho 


(2.3.7) 

(2.3.8) 

(2.3.9) 
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e 


-2R 

o 



e 




o 



Also, the orbital angular velocity. Cl, is given by 


(2.3.10) 


(2,3.11) 


Cl = J Ai/a 3 (2.3.12) 

Now, we can write Eq. (2.3.2) in terms of global unit vectors 
I, J, K as: 






(2.3.13) 


where [T^.] and [T ] are transformation matrices given by Eqs . 
& ot 

(2.2.7) and (2.2.8). Also, 


t = R cos & J + R sin & K (2.3.14) 

o o o 

So, Eq. (2.3.1) becomes, 

& = [R cos & + r_ cos (0 + ot)] J 

Z O fa 

+ [R sin & + r- sin (0 + ot) ] K (2.3.15) 

O L 

where ot is the pitch angle shown in Fig. 2. 

The position vector being thus known, it is easy to calculate 
the velocity vector For, differentiating Eq. (2.3.1) yields 

$ 2 = d$ 2 /dt + dr 2 /dt 

= dt 2 /dt + 2 x r 2 


(2.3.16) 
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where. 


dr. 

4 

dt" 


I 

& + <*. 


K 


r_ cos (<S>+«) r ? sin ( e+ot) 


= - [r, (0+<x) sin ( d+ot )l J 

4 * 

+ [r- (6+ot) cos (£+<*)] K 

it 

And differentiating Eq. (2.3.14), 

dt Q /dt = ( R q cos 6 - R c & sin 6 ] J 

+ [ R sin & ~ R_ ^ cos & ) K 


(2.3.17) 


(2.3.18) 


Therefore , 


$ = [ R cos 

v 2 o 

+ [ R sin 
1 o 


© - R & sin & - r 7 (£+«) sin (©+«) ] J 
o 4 

d - R e cos S + r, (<3>+«) COS (©+«)] K 
o 4 

(2.3.19) 


The .a e nltadea of the poaition and velocity vectora. Ej and V 2 


are, therefore, 


r = ( R " + r_ 

k 2 v o Z 


2 + r„ 2 + 2R r 7 cos « ) 
o « 


1/2 


(2.3.20) 


V 2 - R 2 +R 2 ® 2 + r o 2 ^ + ** ^ 

v 2 o o 4 

+ 2 r 7 6 (1 + «’) <* 0 ^ cos « - R c sin «> (2.3.21) 


and similarly. 


v 2 = r 2 + r 2 e 2 + r 2 & 2 (1 + «’)' 
v i o o A 


& (1 + «’) {-R d cos * + R c sin <*> (2.3.22) 


+ 2 fj » ~ ' ' o 
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&• 3. 2 Firubl orbit pArAmotors 

The equations of motion as obtained in Sec. 2.2, viz. 


or 


f 2e sin 

( l+e cos & J ^ 


_ a i v~ + 1 ) + -s — 7 37 cos « sin « = 0 

•e cos © J l+e cos & 


(2.2.13) 

can be solved numerically to yield a and «’ at any & for any set 
of given initial conditions. Then, from last subsection, and 
R 2 can be calculated and therefore specific energy (total energy 
per unit mass of body M 2 ) «, becomes: 


H . v 


2a, 


R 2 2 


(2.3.23) 


leading to 


13 


R V 

* 2 - *,/ 


(2.3.24) 


and. 


R V 2 2 R V 2 

2 2 ' ' " 


■{[ -V-- 1 ) *( -V-K 


R V 

2 *2 If. 221 2^1 

"2 J cos /3 V 


1/2 


where f5r is the angle between vectors S and , that ij 


cos ft 


$ 1 

*2 2 

R 2 V 2 


(2.3.25) 


(2.3.26) 


2.4 AN APPROXIMATE ANALYTICAL SOLUTION 

The equations of motion, given by Eq. (2.2.13) is a nonlinear 
second order differential equation. To get an analytical 
solution, we must be able to linearize it. Although linearization 
will obviously introduce some error, it is still useful in the 
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sense that we can get an Insight into the basic nature of 
problem. 


ot” 


2e sin 6 
1+e cos & J 


(ct’ + 1) + 


3 

1 + e cos & 


cos ot sin ot = 


( 2 . 


ot' 


2e sin 6 ot ’~j + 
1+e cos d J 


3 

1+e cos & 


cos ot sin ot 


2e sin 6 
1+e cos & 

Rearranging, we get, 

(1 + e cos & ) ot” - 2 e sin & ot’ + 3 sin ot cos ot = 2 e 

(2- 

Ue put, 

(1 + e cos ot = y C2. 

therefore , 

y ’ = (1 + e cos & ) ot’ - e sin & o* C2. 

and 

y» = (1 + e cos 0) ot” - 2e sin & ot’ - e cos ot (2. 

Thus, Eq. (2.4.1) becomes, 

y" + e cos 0 a + 3 sin ot cos ot = 2e sin 6 (2. 

Assuming ot to be small, so that 
sin ot me ot and cos a 4 1 
this equation becomes, 

y” + (3 + e cos & ) ot = 2e sin & (2- 

Substituting ot from Eq. (2.4.2), we get 

~ . (3 + e cos 0)y _ 2e sin & (2- 

+ 1+e cos 6 


the 

0 

2.13) 

sin & 
4-1) 

4.2) 

4.3) 

4.4) 

4.5) 

4.6) 

4.7) 
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For orbits with small eccentricity, e, we can write approximately 
from binomial expansion, 

y” + (3 - 2e cos 0) y = 2e sin 9 (2.4.8) 

which can be further simplified as 

y” + 3y = 2 e sin 9 (2.4.9) 

This is a linear equation which has the following closed-form 
solution: 

y = A cos 4 3 9 + B sin 4 3 9 + e sin 9 (2.4.10) 

where A and B are constants which depend upon initial conditions. 

If ot and a ’ be the initial values for a and a’ respectively, 
o o 

from Eq. (2.4.2), we have 
y (0) = (1 + e) « o 

so that, A becomes 

A = (1 + e) a (2.4.11) 

o 

Similarly, we can find 

(1 + e) a ’ - e 

B = 7T — 2 (2.4.12) 

Putting these in Eq. (2.4.10), we get 

a = ^ — f A cos -4 3 9 + B sin 4 3 9 + e sin 9 J 

1 + e cos 9 l J 

(2.4.13) 


Similarly, 


a ’ = — f - 43 A sin 43 d +43 B cos 43 d 

1 + e cos 9 l 

+ e cos 9 + e sin 0 a j (2.4.14) 

These analytical expressions readily yield the Instantaneous 
position and velocity vectors for the payload as given by Eq. 
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(2.3.15) and (2.3.19). 

The analytical expression thus obtained for the total payload 
energy when evaluated at small intervals around an orbit enables 
us to determine the optimum "cut-off" location ..as. _ well. as. ^ita 
corresponding major axis gain. 

2. 5 NUMERICAL SOLUTION 
2. B. 1 Methodology 

Ue solve numerically the second order differential equation 
(Eq. 2.2.13) using the NAG subroutine D02BHF . The values thus 
obtained are utilized to obtain the position and velocity of the 
payload to be released and hence the major axis of its final orbit 
in release. In addition, we want to evaluate exactly when the 
payload should be released so that the payload major axis is a 
maximum. It may be noted that the determination of the time of 
release comes down to the determination of the position of the 
centre of mass (its true anomaly, &) of the entire tether system 
on the Initial orbit when the payload is released. One such 
method has been described by Amier^^, but it is applied to 
essentially circular starting orbits. Ue now extend it for 
elliptic starting orbits. 

2. S. 2 Numerical Integration 

To solve Eq. (2.2.13) numerically. It was reduced to two 
equivalent first order differential equations by letting 
Y( 1 ) = « 

= a’ 


and 


Y(2) 


(2.5.1) 
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to that we have 


Y’(l) = Y(2) 

Y ’ ( 2 ) = L e - 8in ~ [Y(2) + 1] 

1 + e cos & 

+ T7 s cos { Y( 1 ) } sin { Y( 1 ) } (2.5.2) 

1+e cos 6 

These equations were solved using Runge-Kutta-Mersion method to 
give ot and cn’ for 6 varying from 0 to 360 degrees. Runge-Kutta 
method is quite accurate for non-stiff equations and stiffness 
test shows that these equations are non-stiff so that a reasonable 
tolerance of ~ 10~* and 10 5 was assumed. With <*■ and ot’ thus 
obtained, the position and velocity R 2 and V 2 can be readily 
calculated as given in Sec. 2.3. Therefore, the specific energy 
becomes 



(2.5.3) 


To maximize it, we must set. 


d<£ 

66 


0 


To calculate d<e/d0, we proceed as, 


Therefore, 



d« 

66 


dV, 


2 60 



(2.5.3) 


(2.5.4) 


Eq. (2.3.20) gives the expression for R 2 as, 


R 2 2 = R o + r 2 + 2R o r 2 COS “ 
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L f f erent lat ing yields. 


2 d W 


2R R J& + r, + 2r_ (R /6») cos « 

O O 4 4 O 


+ 2 R r„ (- sin « ) a’ 
o 2 


(2.5.5) 


here dot denotes differentiation w.r.t. time and prime denotes 
i f f erent iat ion w.r.t. 6. Thus the 2nd term on the right of Eq. 


2.5.4) becomes 


dR- 
/j 2 


. r/ i 

R R + * — + r- R cos « - R r„ sin « «’ 

o o 2 2 o o 2 


(2.5.6) 

R , R_, & and their derivatives are obtained from Eq . (2.3.5) to 

O 4 

(2.3.11). Again, Eq. (2.3.21) gives, 

V 2 = R 2 + R 2 & Z + r 2 9 2 (l+«’) 2 

z z o z 


+ 2r„ d (l+o»’) [R 6 cos « - R sin ot ] 
2 x o o 


(2.5.7) 


differentiation gives 2V £ on the left which is twice the first 
term on the right of Eq. (2.5.4). So, we divide each term on the 
right of (2.5.7) by two after differentiation and that gives us 
following five terms: 


F(l) 


R k /& 
o o 


(2.5.8a) 


F < 2 > = + R o & 


(2.5.8b) 


F(3) = r. 


(!+«’) ( 1 + a ' ) + & 2 


(2.5.8c) 
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F(4) 


r 2 ( 2R o k C 1+ot ’) cos ot + R q 9 2 a” cob a 

- 2 

- 9 (1+ot’) ot’ sin ot R + 9 (l+«') cos ot R 


(2.5. 8d ) 


F(5) = - r 2 £ (9/9) R q (l+«’) sin ot + k o (1+ot’) sin ot 

+ 9 R o ot” sin ot + <9 R q (1+ot') ot’ cos otj (2.5.8e) 

Here also, the derivatives of R and 9 are obtained from Eq . 

o ^ 

(2.3.5) to (2.3.11) and thus. 


dV 2 

v 2 = FC1) + F(2) + FC3) + FC4) + F < 5 ) (2.5.9) 

combining Eq. (2.5.6) and (2.5.9), we get the expression for de/d & 
in terms of ot, ot’ and 9. Numerical solution of (2.5.2) gives ot 
and ot’ for various 6. Our task is to find that 9 for which de/dd 
(a function of ot, ot’ and 9 ) becomes zero. Actually, this 
condition alone is not sufficient to guarantee the maximum. It 
gives several local extrema (both maxima and minima) but the total 
number of the local extrema is rather small (as can be seen from 
the next section). Ue numerically obtain the global maximum. 

At 9 = 0, ot and ot’ are known from initial conditions. Ue 
calculate d«/dS using these values at initial point. Now, 
increasing 9 by a small amount, <5, we obtain the new ot, ot’ from 
the numerical integration of (2.5.2) and check for the sign change 
in d«/dd. If there is no sign change, we keep on increasing 9. 
Thus, we can bracket the range over which a sign change occurs. 
Then, the bracketed interval can be shortened by reducing S. But 
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reducing <5 also reduces the accuracy of a and a’ obtained by 
numerical integration. We have to be judicious in the choice of 

j 

tolerance, <5, so that the error stays within reasonable limits. j 

Je have used more than one value of tolerance to get reasonably 
accurate and constant results. 

The above algorithm is used to set up a computer program to j 

compute the various true-anomalies d^’s where the payload should 
be released to reach the maximum or minimum possible altitudes. 

This enables us to obtain the maximum possible increase in the 
major axis of the payload on release alongwith the "cut-off 
location. 

2.6 RESULTS AND DISCUSSION 

Numerical results are obtained for following parameters: 

H = 10 5 kg, = 500 kg, = 500 kg, 1 (length of tether) = 100 

km. Rp (perigee height) = 440 km. The time (or 0) is measured 
from the perigee passage. The radius of the Earth is assumed to 
be constant and equal to Re - 6378 km whereas p 3.986x10 m /s . 

The tether system is assumed to enter its nominal orbit at 

the perigee with initial rotation and initial velocity ’ 

■ ' ■■ ■ ■ | 

denoting the in-plane pitch-motion. The out-of-plane librations, 

. '■ .. . . : j 

such as roll, has not been considered because of its little effect 
on altitude gain as shown in earlier investigations. 

With two tolerances of <5 = 10" 4 and lo' 5 , the program runs 
are performed for different initial configurations of the tethered 

system. The tabulated results for the following configurations 

' . '' ' ■ . ■ ■ ' ■■ ' ’ . \ ■. ■ ■■ ' ■ ' ■ ■■ ' , ■ ■ . ' . ' ; . ■ ' . ' ' ^ ■ ", ' ' ' 

are presented from Table 2.1 to 2.10. 
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(i) 

e 

= 0; 

% = 0 

; « o ' = 0, 0.05, 0.1, 0.2, 

0.5 

(ii) 

e 

= 0; 

a 9 ~ 
o 

0; « q = 5.8, 18, 45, 60, 90 

degrees 

(iii) 

e 

= 0.01 

; V 

* 0; «o * 0, 5.8, 18, 45, 

60, 90 degrees 

(iv) 

e 

= 0.01 


0; « o ’ = 0.01, 0.05, 0.1, 

0.3 

Cv) 

e 

= 0.05 


0; * ' = 0.02, 0.1, 0.3 


(Vi) 

e 

= 0.05 

, a = 
o 

5.8, a ' = 0.01 
o 





CL = 

O 

45, « ' = 0.01 

o 





CL = 

O 

120, a ’ = 0.3 
o 


Cvii) 

e 

= 0.1; 

CL = 

O 

0, « ' = 0.02, 0.1, 0.3 


(viii) 

e 

= 0.1; 

CL = 

O 

5.8, ot ’ = 0 
’ o 





Cl = 

o 

45, ot ’ = 0.01 
o 





ot = 
o 

90, ot ' = 0.05 
o 





a = 
o 

120, ot ' = 0.1 
o 


(ix) 

e 

= 0.5; 

ot = 

o 

0, ot ’ = 0.01 

’ o 



ot = 0, ot ’ = 0.1 
o o 

a = 10, ot ’ = 0 
o o 

ot = 45, ot ’ = 0.001 
o ’ o 

ot =90, ot • = 0.12 
o o 

Tables 2.2 and 2.3 list the results for a circular orbit, for 

which solutions are available from earlier investigations. As we 

8 10 

find, our results are matching well. Bekey and Amier reported 
that if the release occurs while the payload is stable and along 
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than the perigee by about seven times the length o£ the tether. 
This so called "seven-times rule” for the non-librat ing system has 
been completely verified by our results. For circular orbits have 
been larger ’ at & = 0 for systematically aligned with a local 
vertical leads to greater increase in the payload orbits. 

Other notable features are, the gain in altitude, H^, as well 
as the gain in the semi-major axis Aa, are much larger for initial 
orbits with high eccentricity, e. Again for the same e, the gain 
is more in the case of higher « * with system initially aligned. 
Thus Initial pitch velocity is more important than the initial 
pitch angle. 



CHAPTER III 


TRANSVERSE VIBRATIONS OF THREE-BODY TETHERED SYSTEMS 

3.1 DESCRIPTION OF THE SYSTEM 

The tethered system under consideration consists of three 
bodies having masses m ± , m £ and m 3 as shown in Fig. 3.1. The 
length of the tether connecting the end mass m^ and the middle 
mass m^ and the one connecting the end-mass m 2 and the middle mass 
are denoted by 1 1 and 1 ? respectively. The masses of the tether 
segments are m tl and m t2 respectively. The angles and 

represent the inclinations of the tethers to the local vertical. 
The orbital coordinate system xyz is located at the system center 
of mass, as shown in the Fig. 3.1. CM is assumed to follow a 
circular orbit with orbital rate O. lj and 1 2 are independent of 
each other and the system represents a. three-body tethered 
constellation . 

3.1.1 Equilibrium Configurations 

Assuming fixed length of tethers, it can be shown, by 
suitable analysis of the librational motion of the system, that 
there are four possible equilibrium configurations. It can 

further be shown that only one of these configurations is stable. 
In this configuration, the masses are aligned with local vertical, 
as shown in the Fig. 3.2a. Here, the transverse vibrations of the 
system are considered for this stable equilibrium configuration. 
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3.1.2 Klnemailcfi of the System 

As the motion is in the orbital plane, the displacements are 
represented by the inplane coordinates y^ , z^ ; for the three 
masses; i = 1,2,3. Furthermore, since CM is the centre of mass, 
the following constraint equations hold: 


3 


E 

i=i 


m. 


3 

0 , £ a, z. = 0 

1 = 1 1 


It should be noted that the tether masses are being neglected 
for this analysis. In the Fig. 3.2, we see that 

z ± < 0 

z 2 ’ 2 3 > 0 

also, Z 3 - z 1 — 1 ^ 

z 2 = Z 3 + *2 = Z 1 + 1 1 + 1 2 


so that the 2 nd of the constraint equations become: 

®l z l + m 2 (z t + 1 ± + 1 2 ) + ® 3 < 2 i + 1 i> = 0 


therefore, 

(u^ + » 2 + m 3 ' ) Z 1 + (in 2 + ‘"S'* 1 1 + ®2 1 2 = ° 

(m 2 +m 3 ) l ± + o 2 1 2 
Z 1 ®l +,n 2 +ni 3 

If we define mass ratios in the following way: 

m. 

u = — ; i = 1,2,3 

H i m 

with m = mj + »2 + “3 
then 

Z 1 = ~ ( ^2 + X 2 

Similarly, * 2 = " ^1 l l + ip 2 + ^ 1 2 


(3.1.1) 

(3.1.2a) 

(3.1.2b) 
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and 


- P, 1. 


1 *1 H Z A 2 (3.1.2c) 

These equations can be rewritten in the following matrix - 
vector form: 


* £CJ 


(3.1.3) 


where , 


[C] = 


' ^2 ^3 


- P. 


*1 ^3 


l Y, 


= [C] 4 


1 i a i 


X 2 “2 


(3.1.4) 


3.1.3 Tensions in the Tether 

Let and T£ he the tension in the two tethers. Also, the 
distance of m 2 from the centre of Earth is given by 

r 2 = R c + Z 2 

where, R is the orbital radius of CM. 

' c 

The free body diagram of m 2 and m^ are shown in Fig. 3.2(b) 
and 3.2(c). From these we obtain 

2 ** “2 
®2 ( -“ C 2 ) = " T 2 " — 2 " 
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m 2 N ® 

T 2 = m 2 < ° r 2 ' * 


r- < l ♦ r >‘ 2 


expanding the last term in binomial series and neglecting 
order terms: 


2 P m 2 2 ? 

*2“ r 2 “ —2“ d ~ 2 *£♦..) 
Ro c 


using o> 2 = -tL, we obtain 

R J 
c 


T 2 ~ 3in 2 " 2 2 


Putting z 2 in terms of mass ratios from Eq. (3.1.2b) 
T 2 * 3iap 2 1 U 2 [ lj + (fi 1 + f 3 ) 1 2 ] 

= 3 », 2 lj * i (S ♦ 5 ) 1 ] 


Similarly, from Fig. 3.2(b) 


"3 r 3 ^ * T 2 ‘ T 1 - 


}J m. 


which leads to 

T 1 * T 2 + 3 "3 “ 2 z 3 

Putting the values of T 2 from Eq . (3.1.5) and z^ from Eq. (3 
we obtain: 

= 3mw 2 [£ (m 2 + £ 3 ) l 1 + 1 2 ] 


3.2 MODELLING OF TRANSVERSE VIBRATIONS 

Since for vibrating tether, like vibrating string, we 
neglect its mass. So, we define new mass ratios, taking 
account the tether masses. Thus, 


m, 


m 


i = 1,2,3 


higher 


(3.1.5) 


.1.2c) 

(3.1.6) 

cannot 

into 
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n 


and 


tl 


tl 

in 


in 


t2 


1 2 
SI 


(3.2.1) 


where, m = m^ + m 2 + m^ + = total mass 

Now, from Fig. 3.3, assuming small motions, we can write: 


tl 


Y 1 ~ (Yj - y 3 ) + u ! coa a 1 


* y x Cl - <r x ) + y 3 ?! + Uj_ tt 1 , t) 


where, is a nondimensional distance along the nominal tether 

line measured from the mass ra^Fig. 3.3. Hence, 

= 0 at m. and ? , = 1 at m_ 

1 11 3 

Also, Uj = 0 at ? ^ = 0 and = 1 

Similarly, 


y 


t2 


y 2 (1 - <r 2 ) * y 3 « 2 + “ 2 


where, ? 2 is measured from mass m 2 and 0 - “ 1- 

Also, U 2 (0, t) = U 2 (1, t) = 0 

Now, from the definition of CM, we have 
3 2 1 

E ■!>)*! I *ti (K ti d !f i ) = 0 
i=l i=l 0 

which leads to 


<*1 + 1 ®tl 5 y 


+ C« 2 + j « t2 ) *2 + C ®3 + 1 ®tl + 1 ^2^ y 3 


+ m tl U 1 + m t2 U 2 = 0 


where 


U 1 = J u i <**1 and °2 = J U 2 ^ 


Ue now define, 

y l = y l + y 3 ; y 2 = y 2 + y 3 



rhen we can show 


*1 - { "2 + I "t2 } 


■ "tl U 1 - "t2 U 2 


(3.2.2a) 


y 2 = 


( 1 - { "2 + 7 "t2 }] *2 - { "l + 7 "tl } y l 


" "tl U 1 - "t2 U 2 


(3.2.2b) 


y 3 = 


- { "1 + 7 "tl } y l - { "2 * 7 "t2 } *i 


- "tl U 1 - "t2 U 2 


(3.2.2c) 


tl 


= [ 1 - { "1 + 7 "tl } -<l) y l - { "2 * 7 "t2 } 


^tl U 1 * P t2 U 2 + U 1 


(3. 2. 2d) 


y t2 = ( 1 ‘ { "2 + 7 "t2 } ^ 2 ) y 2 * { "l + 7 "tl } y l 


- "tl U 1 - "t2 U 2 + u 2 


(3 . 2 . 2e) 


3.3 KINETIC ENERGY OF THE SYSTEM 

The kinetic energy (KE) is given by 
. 3 


K = 


j E (y, + *i)‘ 

1 i=i 1 1 1 


+ J E J c« ti d (y ti 2 + * ti 2 > 
i=1 0 

For small transverse motion considered here, the KE 
associated with displacement components in the z-direction is 
negligible. However, the potential energy due to this cannot be 
neglected. So, we take into account the contribution due to 
motion in the z-dlrectlon in the calculation of PE in next 
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section. Thus, KE is 

3 2^2 

* “ 7 £ “i ^i 2 + 1 ,E 1 ”ti 1 y tl d <1 


= 7 E t-i V * T A "tl 1 4 ?I 

i=l 1=1 o 

Substituting the expressions for and y ti obtained earlier, we 
get K in the following form: 

1 


m 


b 


*H + *22 y 2 2 * *33 “l 2 * *44 °2 2 * *55 J °l‘ d ' 


+ a 


66 


J u/ d K 2 + 2a 12 y x y 2 + 2a 13 y ± 


U, + 2a 14 y x U 2 


X 

2 V 2 (a 15 6 1 +S 15 J«l i l d «l) +2 


*16 T 1 U 2 


* 2 *23 »2 U 1 * 2 *24 y 2 °2 + 2 *25 y 2 °1 


+ 2 y 2 ( *26 *2 + *26 J *2 u 2 d '2 J + 2 *34 U 1 U 2 

0 


+ 2 a 35 Uj + 2 a 36 », »j * 2 “ 45 - 2 


U„ U„ + 2 a 46 U 2 I (3.3.1) 


] 


where, a^ etc. are 


functions of mass ratios. For example, 


‘11 


Jij ( 1 - 0*1 + 7 "tl 5 ) + "2 * 2 Mtl) 

* (Mj * 7 *tl )2 * "tl j [ 1 ' C "1 + * "tl’-'l] 


+ ^t2 (p l + 2 p tl 5 
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= Oij - 

\ ^tl )2 CfJ l * *2 * *3 + fJ t2 5 + **1 ( 1 ' 2C "l + 

\ ''ti 5 ) 

+ ^tl 

( r (''l * T + + T "tl )2 ) 


- + 

\ ^tl )2 (fl l * "2 + "3 + "tl + ''t* 5 + (,J 1 + T 

*«> 

- 

+ 1 "tl* (2fl l * "V 


= + 

\ "tl )2 * <"l + 1 "tl> ' 2 ( *1 * T "tl )2 


Therefore, 



P 

II 

1 2 

CHJ. * 1 "tl> ^1 + 1 

(3.3.2a) 

Similarly, 


i 

a 22 

^2 + I * t 2> -^2*1 M t2 )2 

(3.3.2b) 

a 33 

2 

*tl 

(3.3.2c) 

a 44 

2 

"t2 

(3. 3. 2d) 

a 55 

**tl 

( 3 . 3 . 2 e) 

a 66 = 

"t2 

(3 . 3 . 2£ ) 

And , 



“ a 12 

( 1 " C ^1 + I (P 2 + 1 AJ t2 ) 

+ #J 2 ( 1 " (Ai 2 + ^ P t2 ) ] (M 1 + 1 ^tl ) 

" ^3 (Ai l + 1 C#i 2 + 2 * J t2 ) 



1 



+ *tl ( ^2 + 1 ^t2> J [ 1 " (AJ 1 + ? "tl> - f 

0 

i ] «i 



3ti 


+ * J t2 + I j [ 1 " C#J 2 + ? * i t2 ) ' *2 ) ^2 

0 

= ^1 ^2 + I ^tZ 3 + M Z c#i l + I ^tl 5 

+ I M tl (P 2 + I ^tZ"* + 2 P t2 C#J 1 + "2 M tl* 

+ C^i + \ ^ tl ) C^ 2 + \ »t2 ) E_ ^1 “ ^2 “ ^3 " P tl " ^tZ 3 
= 2 (^ 1 + | M tl ) C^ 2 + ^ ^ t2 ) " 0*1 + I ^tl 3 C#J 2 + 1 ^tZ 5 

Therefore, 

a 12 * - <"l 4 1 "tl' C>i 2 + T "t2 ) C3.3.2 e ) 

Again, 

a 13 = " *1 [ X " C ^1 + * ^tl 3 ] P tl + ^2 (#J 1 + 1 ^tl 5 *tl 

+ p 3 C ^1 + 1 ^tl 3 ^tl 

- ^ti 2 J ( 1 " + k _ ? 1 ] *1 

0 

+ ^tZ ( ^1 + I ^tl 3 M tl 

= <**1 * T "tl> '•'tl 1 * '! * “3 ‘ "tl + "t2 ] 

2 1 

" ^1 M tl ' ^tl 2 

- (f'i * 1 ♦'tl 5 "tl ' M tl CM 1 + I ''tl 5 = 0 


( 3 . 3 . 2 h) 
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14 

- "l t 1 - (f>l * 1 "tlO "t2 * ^2 CfJ l + I "tl 

* "3 ( "l + \ ^tl 5 ''tZ 

5 **t2 


1 

‘ P tl J [ 1 ~ (#i l + \ ^tl 3 _ K 1 ] ^1 P t2 

0 



+ P t2 (M 1 + 2 ^tl 3 P t2 


= 

( ^1 + \ ^tl 5 ^t2 [ "l + #J 2 + ^3 + ^tl + P t2 3 



^1 2 ” ^tl 2 ^t2 


= 

<*! + I ^tl 5 2 ■ ^t2 (M 1 + \ ^tl 5 = ° 

(3.3.21) 

*15 = 

P tl ( 1 “ (H 1 + I ^tl 3 ) 

(3.3.2 j) 

*15 = 

"^tl 

(3.3.2k) 

a l6 

" *t2 ( ^1 + 1 ^tl 3 

(3.3.21) 


Now, 

a 23 = ^1 ( ^2 + \ **tZ ) ' H 2 [ 1 ' C ^2 + 1 P t2 ) ) ^tl 

+ ^3 ( ^2 + 1 * i t2 ) **tl + ^tl ( ^2 + 1 ^t2 ) P tl 

1 1 

- ^ t2 M tl J ( 1 " ( ^2 + 1 P t2> ' ? 2 J ^2 
0 

- O', + T M t2 ) *tl t Mj ♦ P 2 + "3 * "tl + "tz 1 
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P 2 P tl P t2 2 P tl 



« 0 

(3.3.2m) 

a 2 4 

' "l (fJ 2 + 7 "t2> "t2 ' "2 ( 1 - ( "2 * 7 



* "3 (>J 2 * 7 "t2 ) ^ 1 2 * "tl <" 2 + 7 tJ t2 ) M t2 



1 

- "t2 *12 -f [ 1 ' (fi 2 * 7 "t2> - f 2 ] ^2 

0 



= 0 

(3.3. 2n) 

a 25 

■ - "tl (>J 2 * 7 ^ 

( 3 . 3 . 2o) 

a 26 

’ "t2 ( 1 * (fi 2 * 7 ' J t2 ) ) • 

(3.3. 2p) 

a 26 

= “ M t2 

(3.3. 2q) 

Also , 



a 34 

= "tl *t2 

(3.3. 2r) 

a 35 

2 

= • "tl 

(3.3.2s) 

a 36 

= ~ ^tl ^t2 

(3 . 3 . 2t ) 

a 45 

= " M t2 

(3.3. 2u) 

a 46 

2 

= - ^2 

(3.3.2v) 

Ue can combine feu terms as: 


a 33 

+ 2 a 35 = - P tl 2 

(3.3. 2w) 

a 44 

2 

+ Z *46 = ” **t2 

(3.3. 2x) 
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a 34 + a 36 + a 45 = - "tl ^2 


Putting (3.3.2a) to (3.3.2y) in (3.3.1), we obtain 


( 3 . 3 . 2y ) 


K 

in/ 2 


{ <fJ l + I "tl 5 ' '"t * T ^tl ’ 2 } y 

*{^ 2*1 -^ 2*1 ^ t 2 )2 } V . 


2 (^ + 2 ^2 + 2 ^t 2 ^ y l y 2 


+ "tl J ”l 2 **1 * **t2 X “2 2 «2 

0 0 

- <"tl U 1 + "t2 ° 2> 2 + 2 »1 ( "tl { 1 - <"l + 1 "«>} °1 

1 

- "tl X <1 "l - "t2 ( "l * 1 "tl> °2 ] 


+ 2 y. 


("t2 { 


1 " O', + 


7 "t2>} 


^t2 X *2 U 2 ^2 ” ^t2 ^2 * 2 ^t2 ^ J 


We put. 


n. 


E Qijtt) 
i = l 


and 


(3.3.3) 


u. 


“ E q i2 (t) ^i2 C? 2 ) 


i = l 


i2 2' 


therefore. 


U 1 = J u l d *1 
0 


n l 

= E ‘inC^ J*ii Cf i )d? i 

i=i 0 
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Ue choose, 

^il = Sin * ” ^ 1 and ^i2 = 8in i ” ^2 

bo that, 

, f 0 for even i = 2, 4, 6, ... 


J" -1 2 sin i n ^ d? 
0 


2-j 2 

in 


for odd i = 1, 3, 5, 


therefore. 


U, 


Similarly, 


U. 


Now, 


n. 


E 

i=l , 3 , . 


Hi 

' in 


n. 


i=l , 3 , . . 


n„ 


J u 2 cUf = E ^ii (t) ; 
" 1 1 i=l 


n. 


0 

Again , 
1 


J “2 «2 ■ £ «I2 




J <1 “l = 5 X ( " 1} 42 q il (t) 5 


0 

1 


n. 


J <2 i 2 «2 = £ ( ~ U 42 q “ (t) 


Putting these, we get KE as, 


n„ 


= 7 « [ a U V * a 22 ^ 2 2 + a 12 y l y 2 * i 5 l a i + 2 .i + 2 q “ 
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+ ^ a i+n + 7 t +n +7 q i 7 + ^ E *1,1+2 y i q il 


i “ 1 i+n^+2 , i+n 1 + 2 H i2 


2 i “l i 

+ 2 a l , i+n^+2 y l q i2 + 2 E a 2 , 1+2 y 2 q il 


n l n 2 


+ 2 *2 , i +n . +2 y 2 q i2 + 2 E E a 2+i,n +2+j q il q j2 

x=l ’ 1 i=l j=l 1 


where. 


= <^1 + 3 " (#J 1 + 2 


a 22 “ C/J 2 + 3 ^ i t2 ) “ ( ^2 + 2 M t2' ) 


<*1 + 1 (P 2 + I ^*t2 ^ 


a i+2 , i+2 


p... (l-Sp^/in) ; i = 1,3,5, n. 


x •* 2 y 4 i 6 i • • • » 


a i+n 1 +2 , i+n^+2 


p.. (l-8p. ? /in) ; i = 1,3,5, ..., n. 


i = 2,4,6, .... n. 


1 , i + 2 


1 [ < 1 ~ x + A J t i / 2 ) } ( 2<| 2 / in ) + 42]; i - 1,3, ..n^ 


"ti ( -’ |2) 


* i = 2 i 4 y • • 


a 2 , i + 2 


{ 1 -C^2 + ^t2 /2)} ( 2-4 2 / in ) + -|2]; i = 1,3,..^ 


p t2 (-42) 


; I = 2,4, . - fi. 


a l .i+rij + 2 


-p 2 (^l^tl^ 25 ( 2-J 2/ in ) + 42; i = 1,3,,-n, 


; I ; 2 9 4 , » • fi# 
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2, i+nj+2 


a 


2+i ,n^+2+ j 


f ~^tl °V p t2 /2) C 2 42 /in) + 42 ; i = l,3,..n. 

0 ; i = 2 , 4 , . .n, 

“^tl M t2 ; i and j = 1,3,.. 

0 ; i and j = 2,4,.. 


Hence the mass matrix is given by 

[M] = m [A] (3.3.4) 

where, [A] is a square matrix of size (^+^ + 2) and is given by 


(A) = 


a„ . a 


12 

a l,n 1 +2 

22 

a 2,n x +2 


a 


1 ,n^+n2 + 2 


a 


2 , n^ +^ + 2 


Sym. 


3.4 POTENTIAL ENERGY OF THE SYSTEM 

Potential energy, V, of the system has two parts, namely, 
gravitational and the elastic P.E. 

V = V + V (3.4.1) 

£ * 

Again the gravitational P.E. is given by the sum of the 
gravitational P.E. of three masses and of two tethers. Of course, 
we neglect the tether def ornations . 

3 v m. 2 ^ ® t i 

V A = - E -r — ' £ ~ r— 

i=l |K C + rjl i=l |R C + r.| 

- £ «i !(* c 3 * V < R C 5 + 

1=1 


Ti)] 


• 1/2 
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2 

E 


M “ti [CR c j + r ti> < R c 3 ♦ r t .)3 


- 1/2 


By expanding the right hand side in a binomial series and 
ignoring terms higher than third order, it can be shovn that 


6 


Y~~ ® “ - A " 1 - E “j r. 
c R 2 i=l 1 1 

c 


1 

2 


+ - ^ 


R 


r m. (r . . r. 
i=l 1 1 A 


- 3 (j . r.) ) 


j • E m*. r. . 
E 2 A 11 


2 E 3 
c 


m ti (r ti ' r ti " 3 - r ti^ > 

Uhere the first term corresponds to the orbital P.E. of the 
system while the second and fourth terms vanish by virtue of the 
CM being the local frame origin. For remaining terms writing r^ 
and r ti in terms of y A and z i and after some algebra, we finally 
get : 


V = V + 

g gorb 


2R 


({cmi * j * tl > - o-i ♦ $ ^ti^} " 2 cl - 3cos2 “l 5 

* {^2 + 1 " tZ 5 ' *1 ^ t 2 ^ 2 } " 2 t 1 - 3 "" 2 » 2 > 

+ 2 Cf'i + 1 ('ti’ CfJ 2 + 1 "tZ 5 X 1 1 2 
/ coe Oj - <» 2 ) - 3 »« “j cos “2 } J 
Now, 1 i^l = y l ’ 1 2°*2 = ~ y 2 

so that with a second order approximation for cos ot 1 and cos « 2 , 


i.e. 
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2 

cos m l - 

cos « 2 2 b 1 - « 2 2 
ve have 


2 2 7 * o 

I* (1-3 cos <* ) = - 2 1+3y Z 

1 z (1 " 3 co ® 2 ct 2 ) = ~ 2 x 2 2 + 3 y 2 2 


and 


{ 


COS 


(«1 - 


**2 ^ - 


3 CO 8 ot„ cos 




} 



Utilising these, we get 

2 

V g * V gorb + “5” ({ 4 I "tl> ' '"l 4 I 

4 { *2 + 1 *t 2 > - * 2+1 'W 2 } <- 2 

+ y l 2 { 3 { CfJ l 4 1 "tl 5 ~ *1 * 1 "tl 

♦ t: { ("i * t - *i * i "tP 2 } 

' 4 

+ y 2 2 { 3 { *2 + I fi t 2 ) " *2 * I ^12 

4 ~T^ { *1 + 1 - * 2+1 "t 2 > 2 } 

^ ^ . l ' ' ^ 

- 2yj y 2 \ *i + i ''ti^ *2 4 1 2 ^ 


L h 2 } 

i tl J 2 } (-2 l, 2 ) 

2 ) 

2 } 

2 } 

2 } 
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+ 41 


Now, 


1 X 2 { ClJ l * 7 "tl 5 - ( "2 * 


*u„ ^ 2 


P 1 = 0 


*u. .2 


= \ -T T i ( ] ( *?* ] x i 

rt *4 ■*“ 


*u„ .2 


21 


1 0 


Using the tension expressions given by Eq. (3.1.6) 


+ h AJ t2 ) } 

(3 

in the tethers 

1 2 

2 

J t 2 ( *r 2 ] 
0 £ 

! ds_ 

i ^ 

. <*1 

1 


* — J 1 2 
21 2 0 

J. (3.1.6) 



V e - ( I { < p 2 * *3> 

0 

+ ''l *2 if } [ 

* U 1 ' 
*1 • 

l*** 


1 1 

- , du„ .2 



+ 

J { *1 ^2 T" M 2 (#i l + 
0 

M 3 5 } ( W~ 2 ] 

<*2 

(3.4.3) 

To discretize. 

we use the same relations as in K.E. 

and 

get the 

following: 

n l 

n l 



U 1 

= E ^^t) 

i=i 

= E -l 2 

i=i 

sin 

in ? x 

and 





* U 1 

n . 

* 2 2 
= E <»ii (t) .2 2 

i = l A " 

(1 + cos 2in ? 1 

) 
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therefore , 


f * u i .2 n i 1 

I *7 J ^1 = q il 2(t) TFT J C1 

0 A i ~ 1 in 


+ cos 2 in ) d£ 


but J (1 + cos 2 in df 1 = 1 + ^_ Bin 2in ? I 1 = 1 

0 '0 


n l 

E TFT -*n 2 <« 

1 = 1 i n 


similarly. 


*vu .2 


Jf-rf) « 


n 2 

E TFT ’u 2 <*> 

1 = 1 in 


Hence , 


Siko* f - - - , 'l*'2 1 2 1 1 2 , _ v 

F 1 "l (JJ 2 * "3 5 4 —IT— ) E TFT Oil ct: > 

v 2 y 1=1 1 n 


if 




{ *1 *2 4 p 2 (P 1 * 'V} .Ej TFT q il 2(t: > 


(3.4.4) 


Putting (3.4.2) and (3.4.4) into (3.4.1), we get. 


V 

s 

£ 

IB/A 

y + v = 
e £ 

( b u ^1 2 4 b 



n A 

n 



1 



1 2 

+ i?i b i + 2,i+2 q il 2(t) + b i,n 1 +2,i+n 1 +2 q i2 2 ) 


+ C 


(3.4.5) 
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where C 1. a conat . tera (V^ etc.) which doee net depend 
any of the generalized co-ordinates so they cancel up 
differentiating in Lagrange’s Eq. {[ |^ ] etc. and 


upon 
whi 1 e 


b n = 3 { (fJ i + 1 ^ti’ -^i*i ^ti^ 2 } 
b 22 “ 3 { (fJ 2 + 3 f J t2 ) - <"l + 7 f'tZ 32 } 

b 12 = ” C,J 1 + T ^ 2 + 1 ^t2^ 

r - - - ^* 1^2 1 2 ’j 

b i + 2 , i + 2 3 { ^1 (/J 2 + + — } 

b i+n^+2 7 i+n^+2 = 3 { ^2 + ^2 C ^1 + ^3 5 } 

and [K] = [B] 

where [B] is a square matrix of size n 1 +n 2 + 2 


[B] 



Sym. 



3.5 EQUATIONS OF MOTION 

Using Lagrange's equations: 


d 

dt 



j « 1-2 N 


(3.4.6) 


(3.5.1) 
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where , 


qj : generalized coordinates 
L : Lagrangian = T-U = K.E. - P.E. 

Qj : generalized force (not arising from potential, 
non- conservative) 

For our problem, = 0 


q 


j 



l n 1 ’ q 12’ q 22 


V 




(3.5.2) 


Denoting y 


y? = q. 


and , 
also , 


q li = q i+2 ’ 1 = 1,2 n i 


l 2i ~ q i+n.+2 ’ 1 = 1,2 - - * ft 2 


Ue can write expressions for K.E. and P.E. (obtained earlier) as, 

(3.5.3) 

(3.5.4) 


1 , - .T 

1 {q j 


T = ^ <q,r [H] (q.) 


and V = j ( qj ) T [K] { Qj } 


where [M] and [K] are given by (3.3.4) and (3.4.6) respectively, 
then Eq. (3.5.1) becomes simply, 

[HJ <q) + [K] <q> = {0} (3.5.5) 

which is a standard eigenvalue problem. 

Premultiplying by [M] _1 and putting [D] = [M] 1 [K], we have 

<q) + CD] <q} = (0) (3.5.6) 

The eigenvalues of [D] are simply the squares of the 
frequencies of vibration whereas the eigenvectors describe the 
modeshapes. 
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3.6 RESULTS AND DISCUSSIONS 

Numerical results were obtained for the foil 
conf igurat ion . 


owing system 


Masses 

: ®i 

5S 

10 5 kg 


: ®2 

= 

io 4 kg 


: ®3 

= 

5 x 10 3 kg 

Tether lengths and masses 

:1 1 +1 2 

J= 

10 km 


: ®tl 

= 

10 k * CE\'T 


: ®t2 

= 

50 kg 

Initial perigee height 

: R 

s 

500 k#^- Aft 




a. ix x i* x o- x perigee neigni : K. = 500 ^ ^ M C O i 

For the above values, the mass matrix [If] and the stiffness 
matrix [ K ] were obtained and the corresponding eigen problem, Eq. 
(3.5.6), was solved for two elements each for the upper and lower 
tethers. Thus, six natural frequencies and corresponding eigen 
vectors were obtained. The Figs. 3.4 to 3.10 plot the mode shapes 
for these eigen vectors. In the lowest frequency modes, we can 
see an almost straight tether although the masses have been 


displaced from their original vertical configuration. Since the 
lowest modes dominate in determining the actual shape of a 
vibrating continuous infinite degree of freedom system, we can 
conclude that tether displacements due to transverse vibrations 
are negligible. A few figures also show the plot for higher 
frequency modes. Here, we can see continuous deformtion of the 
tethers. The noticeable feature, here is that the deformations in 


upper and’ lo^rer , t j&thers 


, Tables y3r.l, 

: jl '9 ~ j 7 .. 

various combinations *6f 


are virtually independent of each other, 
list the first few modal frequencies f 
mass ratios. The frequencies appear to 


or 

be 


quite independent of variations in the tether masses, and ®^-2 
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for all practical purposes when a reasonably light tether is 
assumed. The orbital frequency which depends upon the semi-major 
axis of the orbit undergoes relatively minor variations. It is 
interesting to note that the lowest mode frequency Is of the same 
order as that of the orbital frequency. 



CHAPTER IV 


CONCLUDING REMARKS 


4 .. 1 CLOSING REMARKS 

Some important features associated with the two dynamic 
problems involving tethered space transportation systems 
considered here are summarized below: 

4.1.1 Payload Resease 

The first part looks into the effects of the eccentricity on 

the payload raising potential of the tether used for its 

deployment from the space shuttle placed in an elliptic orbit. The 

influence of pitching librations has also been taken into account 

while the roll has been Ignored. Unlike the "forward swing zero 

15 

libration”, approach of Conway that requires trial and error to 
determine the optimal tether cut-off location, a rather direct 
method of maximizing the total subsatellite energy has been 
employed here. The study clearly demonstrates the usefulness of 
placing the dumb-bell tethered space shuttle system in a highly 
elliptic orbit as it enables a substantially enhanced tether 
payload raising capability. In some extreme cases, e.g. , for an 
orbital eccentricity of 0.5, the increase in the subsatellite 
major axis can go upto as high as ~ 50 times the tether-length. 
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Uith some additional pitching disturbances incidental to tether 
deployment or miscellaneous space operations, this figure of merit 
is likely to be even higher. This finding may add further to the 
multitude of spin-offs of the tether concept already on the anvil. 

4.1.2 Tether Vibrations 

Here, the in-plane transverse vibrations of the two tethers 
connecting a three-body constellation in a near-Earth circular 
orbit have been investigated. Attention is focussed only on the 
frequencies of vibrations in the lone useful equilibrium 
configuration aligned with the local vertical. The natural 
frequencies of the first few important modes are obtained using 
finite element approach. The frequencies turn out to be quite 
independent of the tether mass particularly when the tether is a 
relatively light. It is interesting to note that deformation of 
the upper and the lower tether remain virtually unaffected by each 

other . 


4.2 SUGGESTIONS FOR FURTHER WORK 

(1) The study of payload release can be extended by taking the 
deployment of tether into consideration as well as all the 
environment forces into account, such as aerodynamic drag, 
oblateness of the Earth etc. It may also be desirable to treat 
the problem in utmost generality when vibrations and the 

1 ibrat ions both are considered. 

(2) The optimal location in tha payload rel.aae la aearchad only 
during the flrat revolution around the Earth. It la feaalble to 


.. sains by extending this search to a 

further augment the orbit size gai v 

tethered space system. This may 


few more revolutions of the 
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however. Imply a longer waiting time before the "optimal tether 
cut-off” . 

(3) For three body systems, longitudinal vibrations can be 
combined with the transverse vibrations studied presently. The 
resulting equilibrium configuration and stability analysis may be 
useful in the preliminary design of system controllers. The 
dynamic analysis of the system in elliptic orbits appears to be 
yet another logical extension of the present work. 
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Table 2.1 

Optimum Apogee Height and Position of the System 

Typical Plot : e = 0.5, « = 10°, « ’ = 0 

o o 



SI. No. 

& 

r 

(deg. ) 

Aa r 
(km. ) 

H 

m 

(km. ) 


1 

16.19 

2347.14 

18223.33 


2 

136.48 

1199.99 

15149.67 

TOL = 1 x 10~ 4 

3 

170.51 

839.66 

13835.23 


4 

198.89 

900.57 

14099.89 


5 

240.94 

1488.20 

15974.31 


1 

16.19 

2347.14 

18223.34 

TOL = 1 x 10~ 5 

2 

136.48 

1200.00 

15149.85 


3 

170.51 

839.73 

14099.31 
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Table 2.2 

Optimum Apogee Height and Position of the System 

e = 0 . 0 , a =0 
o 




0 

0.05 

0.1 

0.2 

0.5 


o 






& 

r 

(deg. ) 

0 

118.54 

95.44 

72.20 

40.04 

Aa 

(km. ) 

420.00 

535.30 

536.23 

539.42 

565.70 

H 

m 

(km. ) 

740.74 

985.64 

987.40 

994.37 

1043.01 


Table 2.3 

Optimum Apogee Height and Position of the System 


e = 0.0, « Q ' = 0 


(deg.) 

5.8 

18 

45 

60 

(deg.) 

77.80 

42.41 

19.70 

13.40 

Aa (km.) 

531.30 

499.94 

339.52 

222.95 

H (km. ) 

978.10 

917.9 

617-8 

401-02 
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Table 2.4 

Optimum Apogee Height and Position of the System 





e = 0. 

.01, « ’ = 0 
o 



a 

o 

(deg. ) 

0 

5.8 

18 

45 

60 

& 

r 

(deg. ) 

194.48 

76.05 

131.06 

269.72 

238.80 

Aa 

(km. ) 

521.90 

534.50 

508.50 

353.70 

229.30 

H 

m 

(km. ) 

964.30 

1068.70 

1053.00 

764.00 

545.00 




Table 2.5 




Optimum 

Apogee Height and Position of the System 



e = 0.01, = 0 

deg. 


01 * 
o 


0.01 °- 05 

0.1 

0.30 

& 

£> 

(deg. ) 

157.6 114.12 

92.69 

56.81 

Aa 

(km. ) 

522.80 530.05 

535.85 

554.04 

H 

m 

(km. ) 

963.4 1014.4 

1051.24 

1126.70 
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Table 2.6 

Optimum Apogee Height and Position of the System 

e = 0.05, a =0 deg. 
o 


a ’ 
o 

0.02 

0.1 

0.3 

& (deg.) 
r 

115.43 

81.65 

50.61 

Aa (km.) 

507.60 

546.26 

593.30 

H (km.) 
m 

1588.05 

1630.93 

1733.77 


Table 2.7 


Optimum Apogee Height and Position of the System 

e = 0.05 


« (deg.) 
o 


ot 


(deg. ) 


Aa (km. ) 
R (km.) 


5.8 

0.01 


67 .22 
557.38 
1609.80 


45 

0.01 

120 

0.3 

294.44 

357.88 

438.32 

273.10 

1451-26 

1215 . 15 
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Table 2.8 

Optimum Apogee Height and Position of the System 
e = 0.1, a© = 0 (deg.) 


(deg.) 0.00 

Aa (km.) 556.70 

(km.) 2529.29 


342.05 41.70 

720.90 660.80 

2865.90 2673.54 


Table 2.9 

Optimum Apogee Height and Position of the System 

e = 0 . 1 


ot 

o 

(deg.) 

5.8 

45 

90 

120 

O 

» 

0.00 

0.01 

0.05 

0.1 

& 

r 

(deg. ) 

300.00 

17.00 

274.45 

343.87 

L 

(km. ) 

699.28 

430.31 

71.72 

330.50 

H 

m 

(km. ) 

2723. 42 

2296.14 

1609.00 

2125.50 
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Table 2.10 


Optimum Apogee Height and Position of the System 

e = 0.5 


ot 

o 

(deg . ) 

0 

0 

10 

45 

90 

ot ’ 
o 


0.01 

0.1 

0.00 

0.001 

0.1 

B 

r 

(deg. ) 

0.00 

0.00 

16.90 

10.65 

247.40 

Aa 

(km. ) 

2336.22 

2486.97 

2347.14 

1644.73 

163.00 

H 

m 

(km. ) 

18209.00 

18510.70 

• 

18223.30 

16850.95 

13916.00 


Natural Fr equencl ea £ or the Variation o £ Haas 
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Table 3.2 

Natural Frequencies for the Variation of Mass 
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Natural Frequencies for the Variation of Tether Mass 
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